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Summary 

The power output from a normally vibrating strip radiator is 
expressed in alternative general forms, one of these being chosen 
to refine and correct some particular estimates given by Heckl for 
different numberical ratios of strip width to wave length. An 
exact and explicit calculation is effected for sinusoidal velocity 
profiles when the strip width equals an integer number of half 


wave lengths . 



The aspects of radiation from different acoustical sources and moving 
structural elements in particular are a matter of importance to practical ends 
such as control and also to the sharpening of conceptual views on the subject. 
Analytical representations that are effectively suited for wide ranges of the 
characteristic (or scale) parameters may be found only in connection with simple 
and idealized source models; and fully integrated or self-consistent calculations 
which account for the interaction between a source and its surroundings present a 
considerably more difficult undertaking. The purpose of this note is to correct 
and improve some estimates offered by Heckl, in the treatise "Structure-borne 
Sound" (1973), for the radiation associated with bending waves in finite plates. 
Specifically, the problem envisaged by Heckl postulates a normal velocity dis- 
tribution everywhere on a full plane, or a line when the dependence on one of the 
coordinates is suppressed; and admits an explicitly realizable outgoing wave 
solution throughout the half-space, or half-plane, facing the source. The total 
radiated power can thus be exhibited in a quadrature or integral fashion and its 
measure sought in conjunction with assigned velocity profiles. 

Let a time-periodic velocity potential involving two coordinate variables, 

¥(x,y,t) = Re{^(x,y)e ia)t } , 

satisfy the homogeneous linear wave equation in a medium with constant acoustic 
speed c and equilibrium density p ; then a boundary value problem in the half- 
plane -°°<x< co ,y>0 is uniquely posed by the differential equation 
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along with the specifications 
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v(x) , |x| < 2/2 
y=0 0 , | x | > 2/2 


( 2 ) 


and 


ty = f (0) 
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As 


/ 2 , 2 
, r = /x 


x +y “ , x = r cos 0 , y = r sin 0 . (3) 


Definite integral forms of solution to this ummixed Neumann type boundary value 
problem are given by 


1/2 


^(x,y) =4 / (k/(x-x') 2 +y 2 ) v(x') dx' 

- - 1/2 U 
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v(0 dC , y >_ 0 , arg A -£ _> 
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wherein denotes a Hankel function of the first kind and 


v(£) = / e v(x) dx 


(5) 


2/2 
I 
- 2/2 

represents a complex Fourier transform of the velocity or source function. 

Since the time average energy flux vector (* symbolizing complex conjugatior. 

N = - y pu ImtyVijj*} 

is divergenceless throughout the half-plane y > 0 , the total radiated power P 
can be expressed as an integral of the radial component of N over all directions 
in the far field or, equivalently, as an integral of the y component of N alonj 
the source line itself; from the latter follow, on invoking (2) and (4), the al- 


ternative versions 
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1 m 3 * 

P = - -z pto / Ira{i^(x,0) -w- ip (x,0) } dx 

2 -i/2 37 


i/2 

= Re / v(x)H^ (k |x-x* |) v (x' ) dx dx' 
* -i/2 U 

- m / JSslli dc . 


( 6 ) 


Heckl chooses the coordinate function 


/\ . mr , . JL to * & . . i ,,v 

v (x) = v n sm -y (x+-) , n = 1,2,... , v n = v n , - y < x < — (7) 


to simulate the modal velocity profiles of a simply supported plate and substitutes 
the appertaining transform function 


|v(rj l 2 = V 2 
' 1 n 
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( 8 ) 


into the final version of (6) ; he thus arrives at a single integral expression 


2.2 


P - r^C'nv ) / 


n 
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(9) 


for the power output per unit length of the plate parallel to its edges, and 
remarks that "this integral unfortunately cannot be evaluated in closed form; 
therefore, several approximations will be employed here." The latter, differ- 
entiated by magnitudes of the characteristic parameter kf, , include a first 
estimate 

P = i pu v 2 U/mT) 2 , kfc « nit (10) 


which is achieved through a preliminary averaging of the function (8) prior to 
substitution in (6). Another estimate 
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P = pcfc/ , » tlTT 

n 


(ID 


is given, in keeping with appropriate views about the relative importance of 
contributions to the integral (9) near £& ** + nir ; and it is merely asserted, 
as regards the last estimate, that 


1 ? i /? 

P pcfcr (ki/Tr) , for k£, = nir 
b n 


( 12 ) 


with a reference to Lyon and Maidanik (1962). 

If the second of the power representations (6) be adopted, on the other 
hand, and made explicit with the source velocity (7), the consequent double integral 

, ' i/2 

P = -g poav^ / [cos (x-x 1 ) - cos (x+x’+£) ]Jp(k |x-x’ |) dx dx* , 

-■■2V 2 

is reducible to a single integral after effecting the change of variables 


x-x’ x+x’ 

K = , n - 


Ji ft 

and this operation yields a replacement for Heckl’s expression (9), namely 

? 1 . r 
P = pw(v n t) / J Q (k£0 1(1-0* cos nir ^ ~ 3 d£ 


(13) 


or 


k % 


P - ± pcJ^ / J 0 (11)K1 --jjj) cos sinig!] d„ CU> 


involving the zero order Bessel function . The above integrals are preferable 

to (9) insofar as they can be estimated with facility and precision when 

k£ « nir or kS, » nir ; and permit an exact evaluation if k£ = nir , i.e., 

, 2 nTT 

P = 4 pcv^Cnu/k) / J 0 ( y) [ (1 - ~) cos y +~^ sin y] dy 
= || (-l) n (nTrv n ) 2 {J 0 (nTr) - ^ J^nm)} 


( 15 ) 
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where J.^ denotes the first order Bessel function. It is a simple matter to 
recover the estimate (12) by disregarding the second term of (15) and utilizing 
the asymptotic approximation 

Jq(utt) ~ (2/nit)^ 2 cos(n7r - ^) = (-l) n (nir) ^^ 2 , nir » 1 
in the first term. 

Heckl's prediction (11), which is independent of the wave length X = 2ir/k 
and representative of the leading term in a short wave expansion, receives support 
from the analysis of (14); however, the estimate (10) does not adequately dis- 
tinguish between even/odd symmetry of the velocity (source) profile, or parity of 
the integer n , at long wave lengths. To elucidate the latter aspect in mathe- 
matical fashion, a power series development for the Bessel function in (13) is 
indicated, with rapid convergence of the resulting (integrated) terms when k£ « mr 
The measure of power obtained on setting Jq(x) = 1 f° r x « 1 , 

.1 2 ^ 1 

P = — Dtu(v Si) f [(l-£) cos nTi£ + sin nnC] d£ 

h n q nir 

= tt poiCv^J,) 2 • / sin mr£ d£ 

■ = f pw(v Si) 2 1 ... ~. c ?s .. ..nZE , k£ « nTr 

n (nTT ) 2 

possesses an order of magnitude comparable with that in (11) , so long as n is 
odd and the corresponding velocity is symmetric relative to the midpoint of its 
range. If n assumes an even value and the two term approximation Jq(x) “ 1 — 
be employed it turns out that 

1 2 2 

P = -r pw(v Si) (k$,/ntr) , kJt « htt , n even 

o n 
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which implies a significantly lesser output from asymmetrical sources at long 

wave lengths. Thus, the complementary estimates 

2 ■ 
paj(v 1/ niT) , n *= 1,3,... 

P = . n 7 ? k£, « nir (16) 

— pu)(v^£/mO (k£) , • n =» 2, A,... 

* 

supersede the single one proposed by Heckl , and their refinements are readily 
secured through the inclusion of additional terms in a Bessel function polynomial 
approximation. 

The exact formula (15) which manifests the relative predominance of power 
output when the constructive interference condition k£ = nff holds, is directly 
forthcoming once the definite integrals 

x 

/ J 00 sin(x-y) dp = xJ 1 (x) 

0 

and 

X 

/ yJ Q (y) sin(x-y) dy = i x 3 J 1 (x) 

0 

are available. A newer estimate on the basis of (14) , appropriate to the short 
wave lengths specified by the inequality 


e = nTr/ki « 1 , 

rests on those of the functions 


and 


x 

K(x) = / J (y) (x-y) cos ey dy 

A " 


X 

L(x) = f J n (y) sin ey dy 

rv " 


It is apparent from Heckl' s integral (9) that the distinction between the sym- ■ 

metry cases originates in the trigonometric factor sin^ ( - ^A Jl H) } which acquires 

2 2 ^ 
the forms cos (C^/2) and sin (C&/2) according as n is odd or even; and these 

' are evidently of different magnitudes at small values of C& , 
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for x -»■ 00 . Taking account of the representation 


, / \ 2 , sin Cx 

J n (x) - — / — d? 

0 v 

it can be demonstrated that 


K(x) = 


A-z 


+ Re ^ / 
IT 


i r/e ix ^ +E > . ^ ^ 


1 \ (C+e) Z (^-e) 2 y , / /2_ 1 


and 


00 / ix(C+e) ix(C-e) . 

L(x) = Re i / - — ' — ^ 

TT x l C+e C-e 


t 2 -l 


where the principal contributions to the latter integrals stem from the neigh- 
borhood of C = 1 if x » 1 . The deductions made in accord with the latter 
circumstance, viz. 


K(x) * 


and 


L (x) - - 


r'iirx 


1 j cos [x(l+e)-ir/4] | cos[x(l-e)-7:/4] 
/2ttx 1 (1+e) 2 (1-e) 2 


) cos[x(l+e)-TT/4] _ cos [x(1-e)-tt/4] 
1+e 1-e 


x(l-e) » 1 


underlie the estimate 


P k — pc£v 


(_l) n+1 (mr/k£) 2 


n l/7^p /s<k« 3/2 a-@ 2 , 


cos(k£-Tr/4) f , (17) 


kfi, >> nir 


which surpasses the companion one furnished by Heckl. 
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